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We consider a search for the CP-violating angle Sep in long baseline neutrino oscillation experiments. 
We show that the subleading <5cp-dependent terms in the z/ M — > v c oscillation probability can be 
easily obscured by the ambiguity of the leading term which depends on |5m|i|. It is thus necessary 
to determine the value of Sm^ with a sufficient accuracy. The survival events, which can be 
accumulated simultaneously with the v c appearance events, can serve for this purpose owing to its 
large statistics. Therefore, the combined analysis of v a appearance and survival events is crucial 
to provide a restrictive constraint on 5cp. Taking a test experimental setup, we demonstrate in the 
Scp-Sm^i plane that the analysis of v c appearance events leads to less restrictive constraints on the 
value of <5cp due to the ambiguity of Sm^ and that the combined analysis efficiently improves the 
constraints. 

PACS numbers: 11.30.Er, 13.15.+g, 14.60.Pq, 14.60.Lm 



I. INTRODUCTION 

Neutrino oscillation has been strongly suggested by 
an accumulating number of experiments using a variety 
of neutrino sources 0- Q . providing rich information 
on the flavor structure of the lepton sector. The mass- 
and the mixing-parameters of neutrinos are nonetheless 
still not completely known; following the definitions of 
Ref. Q, the unknown parameters include the value of 
the mixing angle #13, the CP-violating angle Sep, and 
the sign of Sm^. The current upper bound on 613 Q 
is expected to be improved down to sin 2 26*13 ~ 0.01 by 
the next generation of nuclear reactor [g, 0, H, |J| and ac- 
celerator experiments 0, 0, 0] , while <5cp and the sign 
of Sm^ are to be investigated by future long baseline 
neutrino oscillation experiments |lCt lll| . 

In this paper, we consider the search for the leptonic 
CP-violation angle by long baseline neutrino oscillation 
experiments using a conventional superbeam. CP vi- 
olation can only be observed through the flavor-changing 
processes such as — ► v c l . It is challenging, however, to 
extract the value of Sep from this oscillation since the u e 
appearance probability is suppressed by the small value 
of sin 2 26*13 ■ The search for Sep is made more difficult 
by backgrounds such as ^-contamination of the incident 
beam, which is inseparable from the signal at the detec- 
tor, and the neutral-current events that are misidentified 
as electron events. Further obstacles to extracting the 
value of Sqp are the uncertainties in the other mixing 
parameters |l2t Il3l \\a\ . Among them, the uncertainty 
in the value of #13, which is bounded only from above, 
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has been widely studied 0, 0, 0] . The value of #13 en- 
ters into the leading term of the v e appearance oscillation 
probability, and the lack of knowledge on its value makes 
it difficult to search for the subleading effect of Sep ■ The 
ambiguity of l^m^l causes a similar problem: the value 
of lawful also enters into the leading term of the appear- 
ance probability and its current experimental uncertainty 
is large enough to obscure the subleading <5cp-dependent 
effect. This uncertainty can also obscure the dependence 
on the sign of Sm^ and make the Sep search more diffi- 
cult. It is hence necessary to constrain the value of {Sm^ \ 
in searching for the CP-violating angle. 

The value of l^m-g-J can be precisely determined using 
survival events accumulated concurrently with the v e 
events 0, 0, 0|. This is due to the large number of 
Vfj, events expected from the large i^-fiux available in 
superbeams. A combined analysis employing survival 
events as well as v c appearance events can thus place 
strong constraints on the values of <5cp and Sm^. We 
carry out an example analysis for a test setup fixing the 
parameters other than <5cp and Sm^. We demonstrate 
that the analysis of the v e appearance events alone does 
not sufficiently constrain the value of Sep in the presence 
of the ambiguity of l^m 2 ^ | and that the combined analysis 
efficiently improves the constraint. 

This paper is organized as follows. In Sec. |nj we use 
an analytic expression of the v appearance probability to 
show that the ambiguity in the leading term can obscure 
the dependence on Sep and on the sign of Sm^ . We point 
out that precise values of |<5™ 2 i| and sin 2 2^13 sin 2 6*23 are 
necessary to search for i5cp • We then use an analytic ex- 
pression for the survival probability to show that it 
can constrain the value of Sm^. In Sec. IIHI we consider 
a test setup and calculate allowed regions of parameters 
in the Scp-Sm^ plane fixing other parameters. We show 
that the combined analysis of v e appearance and sur- 
vival events, in principle, gives a strong constraint on the 
value of Sep and on the sign of Sm^. We conclude our 
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work and give discussions in Sec. IIVI 

II. RELEVANCE OF THE AMBIGUITY OF Sm^ 
IN THE LEPTONIC CP VIOLATION SEARCH 

Assuming that the number of neutrino generations is 
three, we define the mixing angles 0% C {1,2,3}), 

the CP-violating angle Sep, and the quadratic mass 
differences Smh as in Ref. We assume the pa- 

rameter values to be Sm^ — 8 x 10 _5 eV 2 , l^m^l ~ 
(2-3) x 10" 3 eV 2 , sin 2 6»i2 ~ 0.3, sin 2 6> 23 ~ 0.5, and 



sin 2 26*13 ^ (0.1-0.2) 0,0], which accomodate all neu- 
trino oscillation experiments except for the LSND exper- 
iment 

We consider the case where the neutrino energy E and 
the baseline length L satisfies (5m 31 L) / (2E) — 0(1) so 
that the oscillation signal can be observed. We calcu- 
late the — > v a oscillation probability in matter whose 
density p, and hence the electron number density n Cl is 
constant. Taking up to first order in Sm\ x and in a = 
2y/2G ¥ n e E = (7.63 x 10- 5 cV 2 )(/V[gcm- 3 ])(£/[GeV]), 
we obtain 01 
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where Sjj and Cy denotes sin 6*y 



and cos Oij , respectively, and the top and the bottom of the double sign are taken 
when Srrifzi > and 5m 31 < 0, respectively. The probability has an oscillatory dependence on the energy due to the 
sinusoidal factor, while the preceding factor in the first pair of brackets, depending weakly on the energy, gives the 
envelope of the oscillation. The leading terms of the argument and the envelope of the sinusoidal factor are given by 
\5ml x \L/{AE) and 4c 2 3 

s 13 s 23j respectively. The subleading terms of both the argument and the envelope depend on 
<5cp and the sign of 5m 31 . One can obtain the knowledge on the value of <5cp and the sign of Sm 31 through subleading 
effects if the leading terms are known with sufficient accuracy; if not, the effects will be obscured by the ambiguity in 
the leading terms. 

The argument of the sinusoidal factor is determined from the energy which gives the peak of the oscillation probabil- 
ity, or the peak energy in short. We thus evaluate the peak energy, including the corrections from energy dependence 
of the envelope, up to first order in Srn^i and in a to obtain 
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where n = 0, 1, 2, 



a' = V2G F n c = a/(2E) = (1.93 x 10~ 4 km _1 )(/?/[g cm~ 3 ]), and the top of the double sign is 
for 6m 31 > and the bottom for 5m 31 < 0. Inverting the sign of Sm 31 changes the signs of the subleading second 
and third terms in the braces of Eq. (J2J, and the dependence on <5cp appears in the third term. The ambiguity of 
the leading term, which we denote by Ad^m 2 ! |), must be smaller than these subleading terms to determine 5cp and 
the sign of Sm 31 from the observation of -E pe akn- 
parameters and n = 0, by 
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The current experimental uncertainty in |<5to 2 i| is about 6 x 10~ 4 eV 2 2j, which is larger or similar to the above 
critical values. It is therefore necessary to constrain the value of Sm 31 within a smaller interval. 

Similar analysis of the envelope leads to another set of conditions. To determine the sign of Sm 31 , we must impose 
the following condition on the ambiguity of 4cf 3 s 2 3S 2 3 



sin 2 26*13 sin 2 6*23: 
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which is evaluated at E 
is 



^peak,o using typical values of the parameters. To constrain the value of <5cP; the condition 
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The current experimental bound on the value of 6*23 is sin 2^3 > 0.9 nm. The value of 6>i 3 , on the other hand, 
is bound only from above to date. Ignorance of the value of 6*13 is shown to be a major obstacle in determining the 
value of Sep through long baseline experiments 0, 0, Q . A possible approach to this difficulty is to combine results 
from reactor neutrino experiments [a, l3i la> l2| • Future experiments searching for the disappearance of v e from reactors 
are expected to constrain the value of $13 independently of the values of other parameters such as <5cp and matter 
density. Reference || suggests that future experiments can constrain the value of sin 2 2#i3 with < 10% accuracy if 
the value of sin 2 26*13 is as large as 0.1. Such reactor experiments can be performed in advance or concurrently with 
accelerator-based Sep searches. In this prospect, we assume that the value of sin 2 26*13 sin 2 6*23 will be known with 
a reasonable accuracy by the time of 8qy> searches, and we keep the value of sin 2 26^3 sin 2 6*23 fixed in this paper to 
explore the best-case scenario for the long baseline experiments. 

The survival events, which can be accumulated simultaneously with the v c appearance events, can be used 
to constrain the value of Sm 31 owing to the large statistics available. The energy dependence of the survival 
probability is calculated, up to first order in <5to 2 i and in a, as |17| 
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where the top and the bottom of the double sign in Eq. {JJ are taken when 5m 31 > and 5m 31 < 0, respectively. 
Since the value of as 2 3 (l — 2c 2 3 s 2 3 )/(l — c 2 3 s 2 3 ) is negligibly small compared to the leading terms under the current 
experimental limits, the observation of the energy dependence of this mode gives the value of 
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We then obtain two possible values of Sm 31 , one being 
positive and the other being negative; their absolute val- 
ues differ from each other by twice the Srn^i factor in 
Eq. |JSJ|. This constraint on 5m 31 would contribute to 
determining the value of Sqp and the sign of Sm 31 . Con- 
sequently, a combined analysis of the v c appearance and 
survival events is crucial for measuring <5cp- It gives 
the simultaneous constraint on <5cp and Sm 311 and thus 
the result shall be presented in the Scp-Sm 31 plane. 



III. NUMERICAL ANALYSIS OF AN EXAMPLE 
SETUP 

In this section, we consider a test experiment setup and 
obtain constraints on i5cp and Sm 31 from the oscillation 
event spectra, which are numerically calculated without 
employing any approximation formulae of the oscillation 
probabilities. We fix other parameters including 6>i 3 and 
6*23 as mentioned in Sec. [HI We show that the values of 
#cp and Sm 3l are significantly constrained by performing 
the combined analysis of v c appearance and survival 
events. 



A. Example setup and the method of analysis 

We consider the following example setup. A wide band 
beam of neutrinos is produced at the upgraded Alternat- 
ing Gradient Synchrotron (AGS) at Brookhaven National 
Laboratory (BNL). We assume the flux of neutrinos given 
in Fig. ^ which is obtained by fitting the flux presented 
in Fig. 3 of Ref. [18|. Neutrinos are detected by a wa- 
ter Cerenkov detector which has 500 kt of fiducial mass 
and is placed 770 km away from BNL in the Kimballton 
mine, Virginia [l9| . We exclusively consider quasi-elastic 
events vi + n — * l~ + p as signals. We assume that its 
detection efficiency is unity, the data acquisition time is 
5 x 10 7 sec, and the matter density is 2.6 g/cm 3 . 

The expected number of vi events, where I £ {e,/x}, 
within an energy bin Ei < E < Ei + \ is evaluated as 

(N^ - H )) = 

r E i+i fM(E) do-fa) (9) 

TN J e AEe{E)L_^i P{ ^ _> ^)^_, 

where T is the data acquisition time, J\f is the number 
of target nucleons in the fiducial volume of the detector, 
e(E) is the detection efficiency, f( v »>{E) is the incident 
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FIG. 1: The neutrino fluxes used in the example analysis. 
The flux of muon neutrinos (solid line) is obtained by fitting 
the flux presented in Fig. 3 of Ref. [lSj . The flux of electron 
neutrinos (dotted line) is assumed to be 0.7% of that of muon 
neutrinos with the same energy dependence. 



Up flux, and da^/dE is the cross section of the detection 
reaction [23] • We present in Fig. the calculated event 
number spectra of v c appearance events and those of 
survival events. The parameters are taken as I^jti^I = 
2.5 x lCT 3 eV 2 , 8ml x = 8.2 x lCT 5 eV 2 , sin 2 29 12 = 0.84, 
sin 2 26*23 = 1-00, and sin 2 20 13 = 0.06. These figures 
show clear signals of v c appearance and disappearance 
between 1 GeV and 2 GeV. 

It is necessary in the analysis to take backgrounds into 
consideration. Significant sources of backgrounds in the 
v e appearance signal are ^-contamination in the incident 
beam, and misidentification of neutral pions produced 
through neutral-current interaction. In the following, we 
take into account only the f -contamination which is not 
separable from signal at the detector. The number of the 
neutral-pion background events is difficult to estimate 
theoretically since it depends on the model of the pion- 
production process, details of the detector design, and 
methods of the event selection; we leave the consideration 
of the neutral-pion background to a future work. The 
expected number of i>i events in an energy bin Ei < E < 
E l+1 is then given by (A/f () ) = (N^ -> Vl )) + {N t (v c -> 
vi)) where (Ni(v c — » vi)) is defined as in Eq. with the 
initial replaced by v c . We assume in the following 
that the flux of ^-contamination in the incident beam is 
0.7% of the flux with the same energy dependence, as 
shown in Fig. ^ by a dotted line. 

We obtain the constraint on <5cp and 5m 31 by the 
following procedure. We generate the event number 
spectrum N$ (dcp,$m%i) for various values of 6, 



3CP 



and Sm 31 while keeping the other parameters fixed. 
The number of events in each bin N- Ul \Scp,Sm 31 ) 
includes statistical errors and is distributed around 
(N-"'\Scp,5ml 1 )). We carry out hypothesis testing with 

the null hypothesis: "the parameter values are Sq™ and 
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FIG. 2: Calculated energy spectra of v c appearance events (a, 
b) and Vy, survival events (c) for the setup described in text 
with <5m§i = ±2.5 x 10" 3 eV 2 (+ in (a), - in (b), and both in 

" eV 2 , 
The 



(c)). Other parameters are taken as Sm 2 \ = 8.2 x 10 
sin 2 26>i2 = 0.84, sin 2 2f? 23 = 1-00, and sin 2 20 i3 = 0.06. 



Vfj. survival event spectra in (c) are presented only for Sep = 
7r/2 since they depend little on Sep- The error bars include 
only statistical errors. 



is of statistical significance, which we examine by the x 2 
test using 

X 2 M (Sep, Sm 2 31 ) = 



E 

i=l 



N^\ScP,Sm 2 31 ) 



(10) 



where nt>in is the number of the energy bins. The allowed 
region is obtained as a region in which the null hypothesis 
cannot be rejected for a certain confidence level. 



B. Constraints on the values of Sep and Srn^j 



We carry out the analysis illustrated in the previous 
subsection for three cases: using v c appearance events 
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only, using survival events only, and using both u e 
appearance and survival events. In the following, the 
values of the parameters except Sm 31 and Sep are fixed 
to those given in Fig. [21 and we take |<5m 2 1 i truo ' ) | = 2.5 x 

l(T 3 eV 2 and 6$™^ = -tt/2, 0, tt/2, and tt. 

First, we show the analysis of the v c appearance spec- 
trum only using x 2 ^'(ScpjSm^). We present in Figs.El 
and fallowed regions for 68.3% and 95% confidence lev- 
els obtained from the v e appearance spectrum. We see 
that v e appearance events cannot always determine the 
sign of Sm^i and the value of Sep- While there are indeed 
cases in which the wrong sign of Sm^ does not give any 
allowed region in the presented figures (e.g. Fig. yfa)), 
there are cases where the wrong sign of Sm 31 also gives 
allowed regions (e.g. Fig. |U[b)). The allowed region 
in such cases, however, is not a single extended region 
but a group of small isolated spots, which indicates that 
they are due to statistical errors and background. We 
also note that the allowed interval of <5cp is enlarged by 
the ambiguity of Sm^ in some cases such as Fig. |3] (a); 
in this case, the allowed intervals of the parameters are 
S C p ^ (-0.5 ± 0.25) tt and Sm\ x ~ (2.5 ± 0.2) x lO" 3 eV 2 
for the 95% confidence level. 

The allowed regions in Figs. |3] and 0] can be under- 
stood in terms of analytic expressions of Eqs. |QJ and @ 
as follows. We expect that the two oscillation spectra re- 
semble each other and give a small value of x 2 ^^ when 
their peak energies and peak oscillation probabilities are 
equal, i.e. 

Speak = -EpeX ) i U ) 

and 

P(vn -> v e \ Sep, Smli ; E peak ) 

r>( r(truc) r 2 (true) ,-,(true)\ C^) 

where £p ca k an d are given by Eq. (J2J with n = 0. 

Solutions to these equations are shown in Figs. and 
21 by a dotted line for Eq. and a broken line for 
Eq. iJEJ. The allowed values of Sep and 5m\i are well 
constrained around the region where the dotted and the 
broken lines intersect or get close to each other, even 
for cases with the wrong sign of Sm^. It shows that it 
is difficult to distinguish two spectra satisfying the peak- 
matching conditions Eqs. Ijlll) and (|12fl . Difference of the 
two spectra away from the peak energy is not significant 
under the limited statistics and background. 

Second, we present allowed regions obtained from 
survival spectra in Fig. We show the result only for 
<^cp = since the results depend little on <5cp U . 
The allowed region separates into two horizontal bands, 
one with 5m 31 > and the other with 5m 31 < 0. The 
error of the two values of Sm 31 shrinks down to about 
±3 x 10~ 5 eV 2 at the 95% confidence level. This small 
error is due to the large statistics available. In terms of 
Eq. 0, the two bands correspond to the two possible 
values of Sm 31 that give the same value of Eq. © and, 



consequently, the same energy dependence of the sur- 
vival probability given in Eq. JJJ. These two values are 
indicated in Fig. by dotted lines which lie inside the 
allowed regions. 

Finally, we present the combined analysis by evalu- 
ating x 2 — x 2 ^ c ' > + x 2i ~ v ^' '• The allowed regions from 
the combined analysis are presented in Figs. |H1 and \7\ 
It is shown that the isolated allowed regions are elim- 
inated in these figures. The sign of Sm 31 is hence 
well determined, and its absolute value is restricted to 
(2.5 ± 0.03) x 1(T 3 eV 2 , whose small error also reduces 
the error of <5cp down to about ±0.15 tt or less around 
the true value. The synergy of v e appearance events and 
survival events excludes the fake allowed regions with 
the wrong sign of Sm 31 and improves the precision of the 
parameters. The observation of both event types and the 
combined analysis are thus crucial. 



IV. CONCLUSION AND DISCUSSIONS 

We considered a search for the CP-violating angle Sep 
in long baseline neutrino oscillation experiments. We 
pointed out that it is necessary to take the ambiguities 
in l^ml-J, #i3 and 023 into consideration in CP-violation 
searches. Their ambiguities can obscure the dependence 
of the — » v c oscillation probability on Sqp and on 
the sign of 5m 31 . We then showed that the survival 
events can be employed to precisely determine the value 
of 6m 31 , and be combined with v e appearance events to 
improve constraints on the value of Sep and to determine 
the sign of Sm 31 . 

We numerically verified the significance of the com- 
bined analysis for an example setup. We assumed that 
the neutrino beam generated by the upgraded AGS beam 
at BNL is observed for 5 x 10 7 sec by a 500 kt water 
Cerenkov detector, placed 770 km away at Kimballton 
mine. We took into account the ambiguity of l^mfjj and 
fixed parameters other than Sm 31 and <5cp- We consid- 
ered only the ^-contamination in the incident beam 
as a source of background. We obtained allowed regions 
in the Sep-Sm^ plane through a x 2 analysis. The anal- 
ysis of v c appearance events alone led to less restrictive 
constraints on the value of Sep and the sign of Sm 31 due 
to the ambiguity of l^m^-J. In contrast, the combined 
analysis of v c appearance and survival events was ca- 
pable of constraining the values of Sep and l^m^l with 
the errors of ±0.15 tt or less and ±3 x 10~ 5 eV 2 , respec- 
tively, and also of determining the sign of Sm 31 . 

The values of #13 and 623, which we have kept fixed in 
the present analysis, may not be known with sufficient 
precision by the time of Sep searches. The ambiguity 
in the value of sin 2 2#i3 sin 2 #23 can actually be a hurdle 
in searching for the value of <5cp and the sign of (5m|i, 
as discussed in Sec. ITT1 and in other studies 0E1>EI- 
This, however, does not change our major conclusions, 
i.e.: the ambiguity in the value of I^to^I is also a sig- 
nificant obstacle to the <5cp search, and that ambiguity 
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FIG. 3: Allowed regions for 68.3% (gray solid line) and 95% (black solid line) confidence levels obtained from v c appearance 
events between 0.7 GeV and 3.1 GeV. The width of energy bins is taken as 0.2 GeV. Incident ^-contamination and statistical 
errors are taken into account. The true parameters (cross) are taken to be 8m^ tluc ' 1 = 2.5 x 10 _3 eV 2 and J^p 110 ' = —it/2, 0, 7r/2, 
and 7r in (a), (b), (c), and (d), respectively, while 5m,2i and sin 2 26 l ;j's are fixed to those given in Fig. [5] The dotted line and 
the broken line show the solutions to Eqs. 11 II and 1121 in text, respectively. 
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FIG. 4: The same as Fig. El but 8m 2 ^ tTUc) = -2.5 x 10" 3 eV 2 . 



can be diminished by employing the survival events. 
First, note that the ambiguities of sin 2 26\3 sin 2 623 and 
of |<5m 2 i| affect two separate features of the — > v e oscil- 
lation probability, namely the amplitude and the phase, 
as can be seen in Eq. {Tjl. Since they have independent 
impacts on the i5cp-search in our approximation, the am- 
biguity of |<5?7i 2 i| will still have significant effects inde- 
pendently of #13 and #23- Second, the phase of the 



survival probability, Eq. JHJ, depends little on #13 and 
#23 for small #13 below the current upper bound. Hence 
the value of l^mfjj will be tightly constrained from the 
disappearance probability regardless of the ambiguity 
in an d #23 ■ Our considerations are also supported by 
Ref. [lj| , although its authors analyzed a counting exper- 
iment using low-energy neutrinos and did not consider an 
experiment measuring the energy of neutrinos as we did. 
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FIG. 5: Allowed regions for 68.3% (gray solid line) and 
95% (black solid line) confidence level obtained from sur- 
vival events between 0.7 GeV and 4.1 GeV. Incident v a - 
contamination as well as statistical errors are taken into ac- 
count. The true parameters (cross) are <$Qp UC ' — tt/2 and (a) 
Smli tIuc) = 2.5 x 10- 3 eV 2 , (b) <5m^ truc) = -2.5 x l(T 3 eV 2 , 
while drri2i and sin 2 2#i/s are fixed to those given in Fig. [5] 
The dotted line is explained in text. 



of the oscillation. The improvement of the uncertainty in 
$13 must await future experiments as discussed in Section 

nn 

The combined analysis presented in this paper can be 
applied to more realistic case studies. For that purpose, 
it is necessary to take account of backgrounds other than 
^-contamination, especially that from single pion pro- 
duction events. A neutral pion becomes a source of back- 
ground when the two photons from its decay are not sepa- 
rately detected due to the limited resolution of the detec- 
tor and thereby misidentified as a u c appearance signal. 
The estimation of this background requires careful treat- 
ments. It is difficult to theoretically estimate the number 
of single-7r°-production events owing to its dependence 
on how strong-interaction processes are modeled. The 
number of misidentified events among them further de- 
pends on the details of the setup and method of experi- 
ments such as the design of the detector and the criteria 
of event selection. Nevertheless, there is a case study 
claiming that the number of background events from the 
neutral pions can be suppressed to be comparable to that 
from the incident z/ e -contamination [lfj . We expect that 
the analysis presented in this paper shall remain effective 
when such excellent suppression of background is possi- 
ble. We leave a realistic analysis with the consideration 
of 7r°-background for a future work. 



Since the ambiguity of I^to^I has been shown to be 
controllable, the ambiguities of #23 and #13 are the re- 
maining obstacles to <5cp searches. Of the two, the value 
of 6*23 can be determined from the amplitude of th e v, L 
disappearance probability (c/. Eq. Q) @, |H Q 
This hardly interferes with the determination of l^m^ | 
since the value of I^ti^ | is determined from the the phase 
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